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The mechanism of synchronization in the random Zaslavsky map is investigated. From the error 
dynamics of two particles, the structure of phase space was analyzed, and a transcritical bifurcation 
between a saddle and a stable fixed point was found. We have verified the structure of on-off 
intermittency in terms of a biased random walk. Furthermore, for the generalized case of the 
ensemble of particles, a modified definition of the size of a snapshot attractor was exploited to 
establish the link with a random walk. As a result, the structure of on-off intermittency in the 
ensemble of particles was explicitly revealed near the transition. 
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I. INTRODUCTION 

Yu, Ott, and Chen (YOG) studied a transition to chaos 
for a random dynamical system with the random Za- 
slavsky (RZ) map which describes the motion of parti- 
cles floating on the surface of a fluid whose flow velocity 
has complicated time dependence |^,||; the state of the 
system is being sampled at discrete times. It was shown 
that variation of a parameter causes a transition from a 
situation, where an initial cloud of particles is eventu- 
ally distributed on a fractal, to a situation, where the 
particles eventually clump at a single point, whose loca- 
tion moves randomly in all time. In both situations the 
random motion persists permanently, so the concept of 
attractor is inappropriate. Hence, after transient times 
had elapsed, they took a snapshot of the particle distri- 
bution on the fluid surface, and called it a snapshot at- 
tractor They observed that the long-time particle 
distribution that evolves from an initial smooth distribu- 
tion exhibits an "extreme form of temporally intermittent 
bursting" on the chaotic side near the transition. 

After YOC's work, a type of intermittent behavior 
known as on-off intermittency has been reported by 
Piatt, Speigel and Tresser [Q. On-off intermittency refers 
to the situation where some dynamical variables exhibit 
two distinct states in their course of time evolution. One 
is the 'off' state, where the variable remain approximately 
a constant, and the other is the 'on' state, where the vari- 
able temporarily burst out of the off state. It has long 
been thought that the intermittent behavior in YOC's 
work belongs to "on-off intermittency" [sj-p^. However, 
there has been some confusion about this because YOC 
did not investigate the RZ map from the perspective of 
on-off intermittency. Indeed in the Letter j|] and sub- 
sequent review paper they analyzed the RZ map 
by studying the largest Lyapunov exponent, the size of 



the snapshot attractor, and the simple one-dimensional 
contraction-expansion random map model to understand 
the intermittent transition behavior. But the essential 
geometrical structure related to the mechanism of on-off 
intermittency is absent in their analysis. Besides Heagy, 
Piatt, and Hammcl (HPH) commented that the snapshot 
attractor can undergo a form of intermittent behavior 
that is similar to on-off intermittency, but the size dis- 
tribution of the snapshot attractor computed by YOC is 
quite different from the laminar phase distribution they 
obtained in their random walk model [ p^ . 

Meanwhile, Yang and Ding (YD) have studied a noise- 
driven uncoupled map lattice as a spatially extended sys- 
tem [^^. As a special case, they considered a map 
lattice where logistic map is located at each site un- 
coupled with its neighbors and each map driven by a 
common random variable, which is regarded as a homo- 
geneous background. In their work, similar transition 
which comes from the instability of the synchronous mo- 
tion of the ensemble was found. However, on-off inter- 
mittency of the size evolution of the snapshot attractor 
was proposed merely on the ground of laminar distribu- 
tion obtained from numerical calculations since YD were 
unable to map the size evolution of the snapshot attrac- 
tor into a random walk. In this respect, their study was 
incomplete. 

Recently synchronization in a pair of nonlinear sys- 
tems subjected to the common noise is revisited §0-[|. 
Cade and Basu showed that this synchronization phe- 
nomena is indeed physical in certain cases, and for syn- 
chronization of those system the randomness is not vital 
|]l9| . Moreover we have shown that, in this kind of syn- 
chronization, the distribution of the random variable is 
vital and there exists on-off intermittency in the bound- 
ary of synchronization region ||^. Also the intermittent 
behavior of the size of snapshot attractor was discussed 
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briefly. Note that synchronization of random dynamical 
systems has same meaning of the chaotic transition in 
YOC's work. This notion of sychronization is analogous 
to the synchronization schemey proposed by Pecora and 
Carrol, if we assume that random variable is realized by a 
system as a drive system and random dynamical systems 
are regarded as a response systems Q . 

In this paper, we explicitly show that the size evolu- 
tion of the snapshot attractor really has the structure of 
on-off intermittency contrary to the HPH's comment. In 
Section II, we briefly recapitulate the properties of the 
RZ map, and in Sec. lU the mechanism of synchroniza- 
tion of two particles will be revealed by using the error 
dynamics. In Sec. IV, we introduce modified definition of 
the size of the snapshot attractor for the case of ensem- 
ble of particles, and explicitly show that this modification 
leads to the structure of on-off intermittency. Finally we 
present summary and brief discussion in Sec. ^ 



which represents fluid instabilities at low Reynolds num- 
bers. When ^„ is absent, Eq. (||) is often called the Za- 
slavsky map plj . If we consider a strip \y\ < Kq where 
Kq > k{l — e~")~^, one iteration maps this strip into a 
narrow band |y| < Ki < Kq. Thus long term behavior is 
confined to \y\ < fc(l — e^")^^. Note that Jacobian deter- 
minent of the map is g{a) = e~". Therefore, the map is 
contracting by the factor g{a) after each iteration. In this 
paper, we will restrict our interest to the case of fc = 0.5 
as studied by YOC. Given a uniform distribution of the 
ensemble of particles on the surface (two-dimension), we 
can find fractal distribution of particles for a < ac — 0.31 
after several iterations. When we increase a close to cicj 
snapshots show almost one-dimensional fractal structure 
(very thin line structure). And for the case of a > Oc, 
the snapshot collapses to a point (zero-dimension), and 
particles move synchronously as if they were a single par- 
ticle. 



II. PROPERTIES OF THE RANDOM 
ZASLAVSKY MAP 

A generic Z3-dimensional random map can be described 
as follow: 

r„+i = F5„(r„), (1) 

where r„ e is column state vector and ^„ is a random 
variable. So the map F^^ is chosen randomly at each 
iteration according to some rule generating As it 
was introduced in , the RZ map describes a particle 
floating on an incompressible fluid. There is constant 
divergence transverse to the surface, and this divergence 
leads to contraction on the surface. And there exists 
vortical flow with complicated time dependence. r„ — 
(x„, y„)^ (the superscript 'T' means transpose operation) 
describes the position of a particle on the surface aX t = 
nT, where T is a constant sampling time. F^ is given 
by 

Xn+i =Xn + f{a)yn mod 2tt 

(2) 

Vn+i = g{a)yn + fcsni(a;„+i -t-^n), 

where /(a) = (1 — e~")/a, g{a) — e~", and k and 
a are control parameters. In specific, a gives rate of 
constant contraction of surface, k stands for the pa- 
rameter of vortical flow, x is angle variable, and y is 
radial variable. On the other hand, random variable 
gives complicated time dependence of vortical flow, 



* The scheme is based on the fact that certain chaotic sys- 
tems possess a self-synchronization property. A Chaotic sys- 
tem is self-synchronizing if it can be decomposed into sub- 
systems: a drive subsystem and a stable response subsystem 
that synchronize when coupled with a common drive signal. 



III. SYNCHRONIZATION OF TWO PARTICLES 

From now on, we will consider particles which satisfy 
Eq. (^. Because ^„ is describing geometry of vortical 
flow aX t — nT and particles are sprinkled on the surface, 
they have different initial conditions but feel common 
random driving To investigate the synchronization of 
particles, we will begin to study the dynamical behavior 
of two particles rather than the motion of whole ensemble 
of particles. 

Let us consider a replica = (a;J^,y^)^ which obeys 
the Eq. (^ as does r„. To show the mechanism of syn- 
chronization, we consider the error dynamics (ED) of this 
system, e„ = (u„, — rj^ — r„, which is the difference 
of two arbitrarily chosen particle's positions. Then, the 
whole system z„ = [xn, yn,Un, Wn)^ can be described by 

i-n+i -Fjjr„), 
©n+i = G^„(e„), 
where G^,, is given by 

Un+i = Un + f{a)vn mod 27r 

(4) 

Vn+i = g{a)vn + 2/ccos (0„ + siu (^^^) 

with [A mod B] = {A + B/2 mod B) - B/2 and = 
Xn+i + Cn- Here, we have applied mod operation ac- 
cording to the translational symmetry of u„ in trigono- 
metric functions in Eq. (^. Note that the ED is driven 
by (pn, which is the sum of Xn+i from F^,^(r„) and ran- 
dom variable Therefore is not affected by the 
evolution of e„, hence transformed system have so-called 
skew- product structure mathematically Because the 
trajectory near hyperplane e = is governed by </)„, we 
will study properties of the ED regarding the random 
dynamical variable </>„ as a system parameter (p which is 
constant during the evolution, before considering whole 
four-dimensional system. 
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A. Error Dynamics of two particles 

We will investigate the local and global bifurcation 
structure by finding the fixed points of and its phase 
portrait in order to reveal the underlying structure of 
the ED. When we impose condition u„+i = u„ and 
Vn+i = Vn on G^, Vnf{a) should be integer multiple 
of 27r. If Vn = 0, 



(5) 



Then, there are two fixed points. 



and 

TT - 20 



Stt - 20 




for < < TT 
for TT < 6 < 2tt. 



(6) 



When Vnf{a) equals ±.2tt, ±Att ■ ■ •, there is no additional 
fixed points near transition point ac |22|. The e„4-i near 



the fixed point e* is written as follows: 
e„+i - e* = J*(0)(e„ - e* 



(7) 



where J* (0) = {dG^/de)\e* is Jacobian matrix evaluated 
at e*. Two eigenvalues A!j_, X*_ and the corresponding 
eigenvectors e-|_ , e_ can be obtained by solving eigenvalue 
problem of the matrix 3*{(f>). Pairs of eigenvalues give 
behavior along eigendirection of e± , and characterize the 
stability of e* From Eq. (Q), Jacobian matrix is 

given by 



3*((h)=( ^ ■^^"^ 

ykcoscj) g (a) + kf (a) cos (j) 

And its eigenvalues are 

A± = + g{a) + kf{a) cos 



(8) 



±v/(l + g{a) + kf{a) cos 0)2 - 4g(«)}. 



(9) 



Note that A![, becomes unity for = 7r/2 or 37r/2 irre- 
spective of a and k. In Fig. |l|(b) the logarithm of |Aj_| 
is plotted as a function of for a = 0.3 and k ~ 0.5. 
We will use logarithms of eigenvalues as a measure of the 
stability of a fixed point. Stability along the eigendirec- 
tion e+ changes at = tt/2 and — 37r/2. For < 7r/2 
or > 37r/2, the fixed point e* becomes a saddle point, 
at which trajectory is attracted along e_ eigendirection 
and repelled along e+ eigendirection, since ln|A^| > 0, 
In I All < and Aj_ is real. For tt/2 < < 37r/2 and 
ln|A!5.| ^ region of is divided into two cases. 

First case is when logarithms of two eigenvalues are neg- 
ative and different. In this region, e* become a stable 
node, which makes trajectories attracted along e+ to the 



fixed point, because In |A1| < In |A!^| < 0. For the other 
case, when ln|A^| = ln|Al| < 0, two eigenvalues are 
complex conjugated pair whose absolute values are less 
than 1. Therefore, e* becomes a stable spiral (stable foci) 
and trajectory will be attracted to e* rotating around the 
fixed point. 

Similarly for the other fixed point e°, the Jacobian 
matrix is 



J°(0) 



1 

—k cos ( 



/(") 

g{a) — kf{a) cos 



(10) 



and its eigenvalues are 



= i{l + .g(a) - fc/(a)cos0± 



V(l + g{a) - kf{a) cos 0)2 - 45(a)}. (11) 

Same as the above case of e*, in Fig. [|(c) ln|A^| IS zero 
at = 7r/2, 37r/2 irrespective of a and k. However, the 
situation of stability is reversed. As shown in Fig. |l|(b) 
and (c), e* becomes stable spiral or node(saddle) and 
e° becomes saddle(stable spiral or node) at 7r/2 < < 
37r/2 (0 < < 7r/2 or 37r/2 < < 27r). As a result 
Fig. |^(a) shows corresponding bifurcation diagram and 
e* and e° exchange the stability through transcritical 
bifurcation at = 7r/2, 37r/2 respectively. Note that e* 
is constant but changes its stability with the variation of 
0. 

For some value of close to ti/2 or 37r/2, an attract- 
ing closed orbit and a basin boundary between the sta- 
ble fixed point and the closed orbit were found in the 
phase portrait. Indeed, we have obtained a homoclinic 
bifurcation which can be found in a dissipative pendulum 
with constant torque [p4[ . In our analysis, however, we 
have recognized that the existence of basin boundary as 
a result of global bifurcation does not affect our present 
results. 



B. Synchronization and On-ofF Intermittency 

Our main interest in this section is the mechanism of 
the synchronization of two particles. Therefore, we will 
return to Eq. (^, where 0„ changes for every iteration, 
and consider the stability of synchronized state based on 
the study of previous section. The synchronized state 
corresponds to e* in G^^. Because e* = (0,0) is a 
fixed point of G^^ , trajectory starting from Zq = (ro, e*) 
belongs to a class of a particular solution labelled by 
(rp, e*) where Fq is an arbitrary chosen vector. This class 
of trajectories constitute 2-dimensional invariant subset 
AA embedded in the whole four-dimensional phase space. 
Because Eq. (^ contains random variable 0„, it will fill 
the whole hyperplane satisfying e = e* and will be 
identical to the hyperplane. We say a subspace is in- 
variant if a trajectory starting in the subspace always 
remains in the same subspace. 
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Let us consider a situation where a trajectory, which is 
initially located outside A^, and z„ happens to be located 
in the vicinity of A4 for some n. Whether z„_|_i will 
be attracted or repelled from M is determined by the 
Jacobian matrix of subsystem G^^ at e* given by Eq. (j^). 
The Eq. (|) can be mapped into a random walk when 
we take logarithm of absolute value on each side of the 
equation as follows: 



In |e. 



■ri+l 



= In |J* 



Inle^ 



(12) 



Note that we have divided e„ into |e„| and unit vector 
e„ = (cos0„,sin0„)^, where 6'„ = tan~^ (vn/un)- Then, 
one can regard ln|e„| as n-th displacement of the ran- 
dom walk, and In |J*(^„) ■ e„| as its n + 1-th step width. 
Asymptotic behavior of the trajectory near M is deter- 
mined by the average of the step width, i.e., the bias of 
random walk. This bias is called transverse Lyapunov ex- 
ponent h±, which measures global stability of A4, defined 
as 



dct)d0 p{(j),e) \n\3*{cj)) -6(61)1 



(13) 



where p((/), 6) is probability distribution of and 9. 

If < 0, i.e., the average of step width is negative, the 
displacement of the random walk will decrease in time, 
and |e„| will approach to zero in proportion to e"'*-^. 
Therefore M is transeversly stable, and the whole tra- 
jectories, which initially start from outside of M.^ will 
be attracted to M. asymptotically. If h±_ > 0, the small 
distance between the trajectory and M will increase ex- 
ponentially and at last will be affected by nonlinear term, 
which will reinject the trajectory to the neighborhood of 
Ai. In this case M is called transeversly unstable. The 
transition from the former case to latter one according to 
the variation of system paramter has been investigated 
by others, and is called "blowout" bifurcation ||25|"|27||. 

Practically we can calculate h± as follows: 



lim 



lim I J* 

n — >-oo 



lim iyin|J*(0,)-e.| 

ri. — >oo Ti < ^ 



(14) 



/i_L vs a. graph is shown in the Fig. g which coincides with 
the previous work |jl|,|| . 

It is also known that when is slightly positive, some 
dynamical variables of the system can exhibit an extreme 
type of temporal intermittent bursting behavior: on-off 
intermittency |p^ , p8| . The essential ingredients of on-off 
intermittency are : (a) a hyperplane should contain in- 
variant subset and (b) trajectory near hyperplane shows 
additive random walk in log domain depending on the 
local stability of hyperplane along the transverse direc- 
tion of Jv[ |0,|2^- Trajectory in large negative value in 
log domain will be shown approximately a constant in 



real scale, and corresponds to the 'off' state. In other 
case positive or small negative value in log domain cor- 
responds to 'on' state. In the pairs of particles, we have 
shown that M is the hyperplane containig invariant sub- 
set, and ln|e„| exhibits random walk near M. according 
to the the stability of J*. Therefore we can conclude that 
there exists on-off intermittency. On the other hand, ac- 
cording to a result based on the study of a random walk, 
it is known that the probability of laminar length(i) is 
proportional to L~^/^e~^/^=, where Lg is the length at 
which the systematic drift due to the bias and the dif- 
fusion spread of random variable without bias become 
comparable [p|j7p^. Fig. |^ shows characteristic scaling 
of on-off intermittency which P{V) scales as L^i . 

In this section, we have shown the mechanism in 
the synchronization of two particles using skew-product 
structure and stability analysis of and verified struc- 
ture of on-off intermittency. 



IV. SYNCHRONIZATION OF ENSEMBLE OF N 
PARTICLES 

In this section, we will study the behavior of ensemble 
of N particles. As described in previous section, on-off 
intermittency can be characterized by a biased random 
walk. To establish the link with a random walk with the 
size evolution of the snapshot attrctor, let us consider 
the relation between the ED and the size of the snap- 
shot attractor. YOG considered the size of the snapshot 
attractor Sn as the r.m.s. (root mean square) value of 
distance from the center of mass of ensemble as follows: 



\ 



1 ^ 



(15) 



where are coordinates of i-th particle (z = 1,2,..., N) 
and ijn = N^^ ^ ■ y^. However, y„ is not so good choice 
to establish the random walk relation between Sn and 
Sn~i- Since j/„ is defined as average of and 5^ can 
be expressed by ?/„ and (j/n)^, it is difficult to obtain 
Sn-i from the definition of Sn in the limit of Sn —>■ 0. 
In order to use the result from the previous section, we 
consider the modified measure Sn which is the average 
length not from the center of the ensemble y„ but from 
an arbitrarily chosen reference particle position y^^ in the 
ensemble. Our modified measure is defined as 



(16) 



which is r.m.s. distance in y direction, where = 
yli — yn- The behavior of Sn and Sn can be seen in Fig. ^ 
for the case TV = 1000, a = 0.3, and k = 0.5. Fig. |(b) 
shows the amplitude are slightly different, but bursting 
and laminar behaviors coincide. Moreover, as shown in 
Fig. 0(a), as Sn approaches to 0, Sn becomes identical 
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to Sn- Therefore Sn may be replaced with S'„ when one 
investigates the size of the snapshot attractor near syn- 
chronization transition point. Since our measure Sji can 
be decomposed into r.m.s. value of the ED variable w^, 
the relation between Sn+i and Sn can be obtained from 
the relation between Vn+i and w„. In the measure of S'„, 
the existence of the hyperplane containing invariant set is 
guaranteed by invariant hyperplane of two particles sys- 
tem, that is the hyperplane stands for the state when all 
the particles move synchronously. 



h^^ j dy p{z,y)ln\z{l~2y)\, (21) 

where p(z,y) is invariant probability density. This leads 
to the same transition point which is calculated by YD. 
One should notice that the random variable is multi- 
plied commonly to every particles in the summation of 
Eq. (|l9|), and this leads to on-off intermittency of Sn- 

B. RZ map 



Uncoupled map lattice with homogeneous 
background 



We analyzed one dimensional simple case of globally 
coupled logistic maps which was studied by YD |l|,|l6[. 



Let us consider the system of N particles. They are lo- 
cated at , . . . , yn , and their dynamics are described by 



Vn 



(17) 



where z„ = a^„-f 5 and S,n is a random variable uniformly 
distributed in the interval (0,1]. Note that this system 
is identical to the model studied by HPH to explain 
on-off intermittency. When a = 1 is fixed and b is varied, 
they reported that the transition from non-synchronous 
state to synchronous one is found at be — 2.82, and when 
b > be, the size evolution of the snapshot attractor ex- 
hibit on-off intermittency. We will explain this transition 
behavior with our measure Sn and the ED. If we consider 
the difference between and such that 



then, 



52 



Vn+l Vn+l 

z„<(i-2?/,';-o. 



2 = 1 
1 ^ 

-^z^«)^(l-2y;-<)2 



(18) 



(19) 



If we consider the limit — > and neglect higher order 
terms X)i(^ri)'^ Si(^n)^ ^^ar the transition region, 
Eq. (|is| ) becomes Sn+i — |z„(l — 2j/^)|5„. Then, we 
have obtained the following relation by taking logarithm 
on both side. 



In5„+i = In |z„(l - 2y'^)\ + InSn- 



(20) 



This relation can be interpreted such that the size evolu- 
tion of the snapshot attractor in log domain is governed 
by the random walk with a step size ln|2;„(l — 2y^)|. 
Similarly to Eq. ( [l^ ) the transverse Lyapunov exponent 



In the case of the RZ map, we consider the case of 
vl~0. From Eq. (|), 

<+i = + 2fccos (0„ + <+i/2) sin «+i/2) 

= gia)Vn + k cos 0„ sin[(cot 9]^ + /(a))w^] 
-2fcsin(/)„sin2 [(cot 6*; + /(a)X/2] 

~ [g{a) + fc cos (cote; + /(a))]< + 0(«)^), 

(22) 

where cot6'; = K/v], and = [cot 6*; + /(«)]<. 

Therefore, 

i 

^ ^ E{3(") + ^cos(/.„[cot0^ + /(a)]}2«)2 



(23) 



Then, Eq. (23) have different result comparing with pre- 
vious case, because the summation in Eq. (^ does not 
have a common driving. As a result, it seems that one 
could not deduce Sn from the right side. However, we 
have recognized the crucial property that the value of 
cot6'Jj for all i becomes identical when ~ 0, which 
will be explained in detail in the following. That prop- 
erties means that the whole particles aligned in a line 
passing through reference particle, because tan^^ is the 
slope of a line connecting position of the reference par- 
ticle to that of the i-th particle. In Sec. Ill, we have 
shown that e* becomes a saddle or a spiral for most (p. 
When e* is saddle, particles are repelled along unstable 
direction. However, since |ln|A;^|| < |ln|Al||, particles 
are attracted fast along e_ and repelled slowly along e+ . 
Therefore, particles are scattered around along unstable 
eigendirection after a few iteration, and become a nar- 
row stripe along a straight line from the reference parti- 
cle. When e* is spiral, we have shown that the particles 
around e* rotate around the fixed point without chang- 
ing the shape of distribution. Fig. |^ shows the behavior 
of particles near e* for constant cj) when e* is saddle or 
spiral. After particles are scattered along a line, the Ja- 
cobian matrix maps the line passing through the position 
of reference particle to another line passing through that 
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of reference particle. Therefore, we could approximate 
cot 6n = cot 6*^ = cot 6*2 = ... = cot for ^ 0. 

This result is consistent with the finding of YOC that 
Lyapunov dimension becomes unity on chaotic side of 
transition. This represents that ensemble of N par- 
ticles is embedded in one-dimensional manifold. Also 
from an evolution of uniformly distributed ensemble, we 
have found that near the fixed point all the particles are 
aligned in a line from the position of reference particle as 
shown in Fig. |^. Therefore, we can write as follow: 

Sn+i = \g{a) + A:cos0„(cot6'„ -I- f{a))\ Sn- (24) 

By taking logarithm on both side, we find 

lnS'„+i = 

In \g{a) + fc cos 0„ (cot 6l„ + f{a))\ + \nSn- (25) 

As a result, Eq. (|2^) shows a random walk process in log 
domain. The transition point of the size of the snapshot 
attractor is determined by transverse Lyapunov exponent 
hs, which is the average step size in Eq. ( |25| ) as follows: 

1 " 

hs = hm — 7 In |.g(a) + fccos(/)„(cot + /(q;))| . (26) 

i=l 

In Fig. 1^, the hg, which is numerically obtained by us- 
ing the original map (|^), shows the same result of hj^ 
which is calculated in case of two particles. In this case 
invariant hyperplane is a subspace satisfying = and 
mJj = for all i. By the hyperplane and the previous ran- 
dom walk relation, existence on-off intermittency of Sn 
is verified. Furthermore, in the plot of laminar distribu- 
tion of Sm the distribution is not affected by the number 
of particles and shows same scaling of the two particles 
case, as depicted in Fig. |^, which shows the consistency 
of our results. 



V. SUMMARY AND DISCUSSION 

In this paper, we have considered the mechanism of 
synchronization in random dynamical system. As a spe- 
cific example, we have investigated the random Zaslavsky 
map in detail. From the ED of two particles, we have 
studied the bifurcation structure of invariant manifold 
and found transcirtical bifurcation between saddle and 
stable node. Our results are consistent with those of 
YOC, but are more explicit in the perspective of on-off 
intermittency structure. In case of ensemble of many 
particles in the random Zaslavsky map, the long standing 
problem to map the size evoultion of the snapshot attrac- 
tor into a random walk is at last resolved by exploiting 
our modified definition of the size of the snapshot attrac- 
tor, which is slightly different from the previously pro- 
posed one. From this success of mapping into a random 
walk, we have obtained qualitative understanding why 
this manifestation of intermittency has the same critical 



exponent, which is found from numerical simulation as 
that of the typical on-off intermittency. We emphasize 
that this mapping becomes possible when the ED vari- 
able are distributed on a line along the unstable direction 
of saddle. For further investigation, it will be interesting 
to study the other systems, that seem to be difficult to 
map into the form of a random walk. 
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FIG. 1. (a) Bifurcation diagram of the ED. In this figure w„ 
is omitted because fixed point is always 0. Horizontal axis is 
(/!> and its value is same as in (b) and (c). In the diagram, solid 
and dashed lines stand for stable and unstable fixed points re- 
spectively, (b) and (c) show logarithm of eigenvalues (A±, Aj-) 
of Jacobian matrix at e* and e°, respectively. Vertical lines 
at (/) = 7r/2,37r/2 stand for transition points of fixed points 
where In I All and In I All become 0. 
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FIG. 2. /i^(sohd line) and fts (circle) vs a when k=0.5. 
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FIG. 3. Laminar distribution for two particles. Horizontal 
axis is length of laminar in logarithmic scale. Vertical axis is 
logarithm of probability of laminar length in arbitrary units. 
Filled circle means probability for a = 0.27 and filled triangle 
for OL — 0.3. Gray straight line have —3/2 slope. The case of 
a — 0.27 shows exponential shoulder which is typically shown 
in on-off intermittency models. 
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FIG. 4. (a) Horizontal and vertical axis are Sn and S„ 
in logarithmic scale, respectively. Both axes are logarithmic 
scale. Solid dots are data from 1000 particles when a — 0.3 
and grayed line is diagonal line. Data close to zero lies on the 
diagonal line and that means S„ and S„ become identical, (b) 
shows time series of Sn and Sn for same period. Although am- 
plitudes are slightly different each other, timings of bursting 
and laminar period coincide. 
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FIG. 5. Evolution of ensemble of 5000 particles for constant 
when a = 0.3 and k = 0.5. Left column of three graphs are 
initial evolution of ensemble with = 0, where fixed point at 
e* is saddle. We start with random initial points distributed 
within -0.01 < mq/tt < 0.01 and -0.01 < uo < 0.01 uni- 
formly in order to show the behavior near saddle. From up- 
per to lower graphs, shape of ensemble stretch along unstable 
manifold as time flows. Contrast to left ones, right column are 
that of ensemble when <jf> = tt and fixed point is spiral. The 
distribution, which is initially located in a ring, rotate around 
origin with contracting its area. Note that in these graphs co- 
ordinates are u^n ^'iid Vn/n, i.e. relative displacement from the 
reference particle. 
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FIG. 6. Snapshot of 5000 particles when a = 0.3 and 
k = 0.5. Initially, particles are distributed in the region 
-0.5 < uo < 0.5 and -0.5 < vo < 0.5. (b),(c) and (d) shows 
a fractal distribution of particles when n = 2080, 2088 and 
2092 respectively, (a) shows small area of snapshot attractor 
for n = 2080, 2081, . . . , 2099 in a same graph. In this region 
particles are aligned in a line for each instant time when par- 
ticle are close to origin. Each line stands for a snapshot for 
a given instant time. Note that in these graphs coordinates 
are Un and Vn, i.e. relative displacement from the reference 
particle. 
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